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Abstract 
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1 Introduction 


Study of duality properties in statistical mechanics and quantum field theory models is im¬ 
portant method for non-perturbative investigation of their phase diagram and held content. 
Duality transformation was discovered by Kramers and Wannier [1] in the two-dimensional 
Ising model. Kadanoff and Ceva [2] generalized the Kramers-Wannier duality relation to the 
nonhomogeneous case (the coupling constants are arbitrary functions of lattice site coordi¬ 
nates) with spherical boundary conditions: 

{l[smh2K^{r))-^/^Z[K] = {l[smh2K,{r))-^^^Z[K], fx = x,y, ( 1 ) 

sinh2Kra;(r) • sinh2K_y{r) = 1, smh.2Ky{r) ■ sinh2K'_a;(f) = 1. (2) 

We denote site coordinates, functions and functionals on the dual lattice by ’’tilda” : r, 
cr(f), Ky_{r), H[K,a], Z[K], .... A site coordinate on the dual lattice coincides with a 
coordinate of the plaquet center on the original lattice: r = r + {x + y)/2 and coupling 
constants K_y{r) = K^{r — V) {V = x, y are the unit vectors along the horisontal X and 
vertical Y axes). 

As was already mentioned in [1,2], relation (|^) can not be understood literally. So, for 
example, using the method of comparing high- and low-temperature expansions for deriving 
duality relation ([^) in the case of the periodical boundary conditions, it is hard to take into 
account and to compare the graphs wrapping up the torus. In fact (0) is correct in the 
thermodynamic limit (for the specihc free energy). However for the nonhomogeneous case 
the procedure of thermodynamic limit is rather ambiguous. In [2] this duality relation was 
obtained for spherical (nonphysical for the lattice) boundary conditions. 

In [3], using global Bianchi identities for link formulation of the lattice spin systems on the 
hypertorus, the contributions of the link variables on the topological nontrivial loops on the 
hypertorus was selected in the partition function, but the duality relations was not formulated 
in obvious form in this case. 

Since duality is a popular method of non-perturbative investigation in quantum held the¬ 
ory and statistical mechanics (for review see [10]), it is important to formulate a duality 
transformation for hnite systems. Recently, we have suggested [4,5] exact duality relations 
for the nonhomogeneous Ising model on a hnite square lattice of size n x m wrapped on the 
torus: 

n(sinh2i?^(f))-'/"z(?'-?-)[ir] = in(sinh2ir^(r))-'/" ^ (3) 

r.u, Px-!Py=^ 
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Here are partition functions of the Ising model with corresponding combinations 

of the periodical {px,Py = 0) and antiperiodical {px,Py = 1) boundary conditions along the 
horizontal X and vertical Y axes: 

Zip^,py)[K] = ^exp(^ir^(r)a(r)V(f‘'V(r)), (4) 

H 

and 


/ 1 1 1 



V 1 -1 -1 IJ 

where r = {x, y) denotes the site coordinates on the square lattice of size nxm, a; = l,...,n 
y = 1,... ,m; a(r) = ±1; Kx{r) and Ky{r) are the coupling constants along corresponding 
axes. The one-step shift operators V^;, act on a{r) in the following way 

Xx(T{r) = a{r + x), Vya{r) = + y). (6) 

For the periodical (antiperiodical) boundary conditions along X and Y axes we have 

X^J^^a{n,y) = {-y^a{l,y), X^y^^a{x,m) = {-Y«a{x,l). (7) 

In Ref. [4] the duality relation (4) was proved for homogeneous and weakly nonhomoge- 
neous distributions of the coupling constants. We also have checked the duality relation (4) 
for lattices of small sizes by direct calculation on the computer. As a corollary of (4), we 
obtained [4,5] the duality relations for the two-point correlation function on the torus, for the 
partition functions of the 2D Ising model with magnetic helds applied to the boundaries and 
the 2D Ising model with free, hxed and mixed boundary conditions. 

In this paper we formulate method of derivation of the duality relations for two-dimensional 
Z(A^)-symmetric spin models on hnite square lattice wrapped on the torus. As example, the 
duality relations for the vector Potts model (the N = 2 case is considered in detail) and 
the Z(A^)-Berezinsky-Villain model [6,7] are obtained. Without taking account of boundary 
conditions duality relations for these models was obtained in [8,9] (for review see [10]). In 
principle suggested method it is not hard to generalize for lattices with larger dimensions 
compactihed on the hypertorus and the lattice models with continious global or gauge sym¬ 
metries. 

To formulate the method let us introduce dehnition of magnetic dislocations connected 
with boundary conditions, ’’topological” charge of dislocation and gauge transformations of 
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coupling constant configurations for the vector Potts model. The hamiltonian of this model 
one can write in the following form: 

- + c.c.) (8) 

r,u 

where a spin variable takes N values: a{r) = exp(i^/(r)), /(r) = 0, N — 1, u = x,y and 

Px = P, Py = Q {p,(l = iV — 1) designate the cyclic boundary conditions for one-step shift 
operators (^): 

V^^V(n, j/) = e*^^cr(l, ?/), m) = 1). (9) 

These conditions have the following form for variable l{r): 

l{n + l,y) = + p, l{x,m + 1) = l{x,l) + q. (10) 

For the periodical boundary conditions we have p = 0 and q = 0. 

Using (P), Hamiltonian one can write as Hamiltonian Ff[iP, d, a] with 

the magnetic dislocation hIP’?) corresponding boundary conditions (p, q) and with periodi¬ 
cal boundary conditions for spin variable a{r): 

^ r,ip 
2'jr 

= H K^{r) cos —{AJ{r) + (11) 

r,z/ 

Here = V® - 1 is difference derivative with the periodical boundary conditions, vector 
helds Kj,{r) and determined on the lattice bonds, it is convenient to consider as the 

module and the phase of the complex coupling constant. The magnetic dislocation is 

determined by the phase 

= (4(r),dy(r)) = (pdi?^(r), gdBy(r)), (12) 

which is unequal zero along the boundary cycle Bx and By, setting the space conhguration 
of the dislocation on the torus: 

^Bx{r)= Y. 6By{r) = Y (13) 

r'^Bx v'^By 

where 5^(r — r') is Kronecker 5-function and 

Bx = {{x,m),x = 1,..., n} , Hy = {(n, p), p = 1,..., m} . 

The phase d(f’'?)(r) one can consider as density of a ’’topological” charge Qi, of the magnetic 
dislocation. This charge, for example, for dislocation is equal 

Qv = ( 14 ) 
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We will call magnetic dislocations {p,q = 0, ...,N — 1) as basic magnetic dislocations. 

Note that periodical boundary conditions {p = q = 0) along all cycles of the torus correspond 
the absence of the magnetic dislocations. Nevertheless, for convenience we have introduce 
denotion for this case. 

Hamiltonian (|^) has invariance relative to Z^v-gauge transformations [11] 

^a(r), K'Jr) = (15) 




a [r] = e 


where 0(r) has the periodical boundary conditions. This invariance gives the following relation 
for partition function: 

Note that the gauge transformation of l{r) and d^{r) in Hamiltonian ( PTj ) has form: 

l\r) = l{r) + 0(r), d'^{r) = + A^0(r). (16) 


These transformations lead to both the deformation of the basic magnetic dislocations and 
the appearance of new closed dislocations. Then d'^{r) describes the held of closed magnetic 
dislocations on the torus. It is obvious that the topological charge does not change at the 
gauge transformation. For example, for Hamiltonian d, a] with dislocation we 

have 

Q', = E4‘”('') = + E 

r r r 

Here the periodical boundary conditions for 0(r) lead to vanishing of the second term and 
Q'fj, = Qim- From here it follows that the set of coupling constant conhgurations {[K,d^^^]} 
(contained closed dislocations) one can divide on the gauge-nonequivalent classes with 
corresponding value of topological charge = {pn,qm). Elements of class one can 

generate with help of gauge transformations ([ISD from the basic magnetic dislocation 

Let us briehy formulate the idea of suggested method. Using the Fourier transformation 
method for derivation of the duality relations, we obtain the expression with ^-functions. The 
solution of the corresponding system of equations dehnes the relation between the initial and 
dual spin variable. Usually, for example, see [10], omitting the problem of taking account of 
boundary conditions, the only one solution of this system of equations is written. However for 
the lattice model on the torus we can End many solutions of this system. These solutions one 
can classify over the gauge-nonequivalent classes of coupling constant conhgurations for 
the dual model and also each class has dehnite value of topological charge = {pn,qm). 
Therefore at dual transformation of the partition function it is necessary to sum over all the 
gauge-nonequivalent classes on the dual lattice with coefficients, depending from boundary 
conditions on the initial lattice. 
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2 Vector Potts model 


Now, using the method discussed in previous section, we derive the duality relation for the 
vector Potts model. Partition function (|TI|) of this model one can represented in the following 
form [9,10]: 


[n [Zl ) 


.271 




[d [t] 


r,/i 


Ott 

^exp(-(3A|i] +i—-p)), 

[t] " 


( 17 ) 

( 18 ) 

( 19 ) 


r,fi 


where 


N-l N-1 

i:=n(i:). i:=n( i: )■ 

[i| ■■ ip).0 111 ■■ t„P).0 

In (|T^) we have made the Fourier transformation to vector held (t^(r) = 0,1, ..., N — 1). 
—f3H[t] is Fourier-transform of Hamiltonian (11): 

M 


277 


-mA = T,T,sriK)cos'‘—t^{T). 


k=0 r,^l 


( 20 ) 


Here M = N/2, if N is even and M = [N — l)/2, if N is odd. In (|^ ^v(’S) is Kronecker 
^AT-function mod N: it is one if s = iVL, where L is integer and zero in another case. 


In order to get rid of (5Af-function in ([1^) it is necessary to solve equation 


- /I) = 0 /. mod N. 


( 21 ) 


Nontrivial solutions of this equation on the torus one can write in the form: 

^^ (r) = {r-V) + (h - z/), (22) 

where index a numerates the solutions, /(f) = 0,1, ..., iV — 1 is dehned on a site of the dual 
lattice and aA> (f) is density of the topological charge (corresponding given solution a) of the 
held of the closed magnetic dislocations on the dual lattice 

II - r'), 4“^ = 0,1, ...,iV - 1. (23) 

i&Za r'eFi 

Here by analogy with (0). (0) <iy(f) is dehned on bonds of the dual lattice. For convenience 
we have written the dislocations on the dual lattice by means of closed paths Fj on the original 
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lattice. Za in (^) denotes subset of the paths (corresponding to solution a) from set F of all 
closed paths on original lattice (Fj G F). Vector a^(r) = e^(r) if the direction of circuit over 
path Fj in site r (the direction of circuit is counterclockwise) coincides with direction of the 
unit vector e^(r) = p in this site, otherwise a^(r) = —e^(r). 

Expression (^) it is not hard to obtain, observing, that solution (p^) satishes by equation 
(Pip on site f when 

- p - F) = 0. 

This equation becomes the identity if m the following conditions are fulhlled: 


-y) = -x-y), -x) = -x-y)- 

d^y\f -y) = d^^\f - x), d^y\f -x-y) = d^^\f -x-y); 

4“^ — y) = —d^x^ if — X — y), 4“^ (r — x) = dSf^ {f — x — y). 

Consistency of these solutions on some set of sites requires that these sites belong to closed 
paths Fj on the torus, that is these solutions must be ’’glued” in order to form the closed 
magnetic dislocatuons. 

Let us denote by coupling constant conhgurations on the dual lattice corresponding to 
the solution (|^) . Depending on the number of the solution these conhgurations cointain both 
the closed dislocations non-enveloping of the cycles of the torus and the dislocations enveloping 
of ones. The dislocations of the hrst type remove by means of gauge transformations (pT]) 
on the dual lattice and the dislocations of the second type can be trasformed to the basic 
magnetic dislocations This means that all conhgurations can be classihed in the 

gauge-nonequivalent classes with topological charge 


f 

where p, g = 0,1, ...,N — 1. 

Since the duality relation connects the partition functions, which are the gauge-invariant 
quantities, at removal ^Tv-functions in (|T^) we must sum over the gauge-nonequivalent solu¬ 
tions of equation (|^): 


- P) + - P), 


( 24 ) 


where is dehned on the dual lattice by relations similar to (0)-(0)- Substituting these 
solutions in (|^, we obtain 

1 

N 


z'r\K,d\ = +v^’i) 


p.i? [i] 
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27r 


exT 


- p) + diP’^\r - 9)]}. 


r,fi 


Here we have introduced factor 1/N as taking into account relation (^), it is not hard to 
note that the sum over configurations [/] in N times more than the sum over [t]. Remarking, 
that 

- z/) = 0, 


r,fi 


relation (P^) one can write in compact form 


7{P>Q) 


.271 




N 


N 


p,g 


r,fi 


^ H exp(i^(pg - qp))Z^^’'^\k, d], 

p,g 


(25) 


where 


d] = X: d']) 

[i] 


is the partition function of the model on the dual lattice. 

Let us in detail consider the case N = 2. Here Hamiltonian coincides with Hamiltonian 
d^) of Ising model. In this case from ( pO]) one gets 

= Y.W°KI<) + ai'-\K) cos7r(Afc!(r))|. 




In order to hnd coefficients g^MR) we use the inverse Fourier transformation 


ex 


^ 97r_ 1 97r 97r 

v(^g^^\K) cos^ —i^) = ^ H exp(R:^cos—n - i—nU). 


N 


N 


N 


(26) 


k ' n=0 

Here /i 7 ^ n. Hence it is not hard to obtain for N = 2: 

_ 1 sinh 2 iFy(r), = tanhiFy(r) = 

where the last relation coinsides with (p]). Using these relations, duality relation (p^) one can 
represented in the form 

1 1 

n(sinh2iF^(r))“^/^Z(^’’P[iF] = - n(sinh XI (27) 

r,P ^ r,^i P,q=0 

where 

Z<?.B|A'] = 5 ;: exp eosx(A5e;») = ^ exp i:(A-,(f)p(f) V;fe)ff(f)), 

[1] r,tM [a] f,tJ. 

and a'(r) = ±1. It is easy to verify that (P7D coinsides with duality relation (pp, since the 
matrix 

Px,Py p,q ' 
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3 The Z(A^)-Berezinsky-Villain model 


Now let us consider the duality relation for the Z(iV)-symmetric Berezinsky-Villain model. 
Partition function of this model one can write in the following form [6,7]: 


zlrdi-Ri-E 

[i] 




EEIlexp 

[i] [fc] UM 





(29) 


where 


7V-1 


oo 


i:=n( sz). i;=n( i; )■ 

[/] /(r)=0 [/c] kfj,{r)=—oo 

Here /(r) = 0,iV — 1 is on a site of the square lattice, index {p,q) dehnes the boundary 
conditions din) ans the sum over guarantes the periodicity of the Hamiltonian relative to 
shifts / —>■ /(r) + NL{r), where L is integer. By analogy with the partition function of the 
vector Potts model (1^) one can rewrite in term of the basic magnetic dislocations 


ZBv\K,d] = ^^nexp{-^77;2(r)£(A^/(r) + djf’‘')(r)) -27rfc^(r)]^|, (30) 

[l] [k] L / JV J 


where /(r) satishes the periodical boundary conditions and the density dfj,{r) of the topological 
charge is determined by relations (12)-(|l^). 

For derivation of the duality relation let ns make the following transformations with (^): 


ZBv\K,d] = + di,{r)) - 27rk^{r)Y \ = (31) 


[l] [k] 


T',k‘ [s] [fc] 




r,/i 


IlE^EE / ^^nexp|-^i^/.(r’)[(A^0(r) + ^d^(r))-27rfc^(r)]Viys(r)0(r)| = 


.N 


.i2vriF(r)^ T.Y. DO^exp 


N 

+i—s{r)e{r) \ = (]d 


r,fi 

N .1 

I 2 


1 Ptt 


EEIlexp 


r,fj. 2vriF^(r) 

n '^(E - h) - Ns{r)), 




(32) 

(33) 


(34) 


where 


For derivation 


E=n( E ). /B«=n 

[t] r,ii t^(^r)=-oo 4- 


r2n 


d9{r) 


0 27r 

we have used the summation formula 


O'TT ^ Ott ^ 

E = E 0 < 0 < 2vr, 


N 


z=o 


N 
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the identity 


^ exp[--ir(/-27rfc)^] = 


k=—oo 

and the dehnition of 5-function 


2 ''^^ 2 ir' 


E exp(-—+ it/) 


r2Tr 




Jo 27r 

In order to take off the 5-function in (^) it is necessary to solve the equation 

y A„i„(r- p) = /Vs(r). 


(35) 


Analysis of the solutions of this equation is similar to analysis of (^lf) - (^4f) and leads to the 
followig expressions for the gauge-nonequivalent solutions 


r -V) - Ne^yKif - V). 

(36) 

■ P-p). 

(37) 


These solutions are the basic magnetic dislocations on the dual lattice in corresponding 

the gauge-nonequivalent classes with the topological charge = (pn, qm) (p, q = 

0,1, ■■■,N — 1). Then, taking off 5-functions in (^), it is necessary to sum over all these 
solutions. In result we obtain the duality relation for the Z(A^)-Berezinsky-Villain model 

(11 ^E - ?p))(n d\, os) 

r,fi 2V I\ p~ I\ ~ I\ 


r,/i 


where 


K^{r)K_^{r) =(^)^ h ^ 


As it is shown in [8] this model at = 2 corresponds to the Ising model, what is consistent 
with our result (^), which coincides with (3) in this case. V.S. thanks Dr. A. Morozov for 

the hospitality and the exellent conditions at ITEP, where this paper has been hnished. 
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